Going Nonlinear with Dark Energy Cosmologies 



Eric V. Linder 

Physics Division, Lawrence Berkeley Laboratory, Berkeley, CA 94720 

Martin White 

Departments of Physics and of Astronomy, University of California, 
and Physics Division, Lawrence Berkeley Laboratory, Berkeley, CA §4720 
(Dated: February 4, 2008) 

We propose an efficient method for generating high accuracy (~ 1%) nonhnear power spectra 
for grids of dark energy cosmologies. Our prescription for matching matter growth automatically 
matches the main features of the cosmic microwave background anisotropy power spectrum, thus 
naturally including "CMB priors" . 



I. INTRODUCTION 

The formation of large scale structure in the universe is 
a central element of astrophysics. Apart from the struc- 
ture itself, it also offers a incisive tool for probing cosmol- 
ogy and particle physics. Since structure formation is a 
competition between expansion and gravity, the growth 
as a function of redshift is intimately tied to the expan- 
sion of the universe. Studies of matter density fluctua- 
tions can thus shed light on the physics responsible for 
the accelerating expansion, i.e. dark energy. The growth 
of large scale structure also provides one of the few tests 
of general relativity on cosmological scales. 

One of the most informative descriptions of large scale 
structure, and the most widely used, is the power spec- 
trum of density fluctuations. The power spectrum can be 
probed through observations including galaxy and clus- 
ter redshift surveys, weak gravitational lensing, baryon 
acoustic oscillations, the Lyman alpha forest, etc. To 
translate the measured power spectrum to a cosmologi- 
cal model, one needs accurate predictions for the power 
spectrum in a variety of dark energy cosmologies. This 
typically comes from large simulations, and fitting forms 
derived from such simulations. However the multitude 
of possible dark energy cosmologies is vast, and the nec- 
essary simulations are computationally expensive. The 
first dynamical dark energy simulations appeared only 
in 2003 1, 2] and even constant dark energy equation 
of state model simulations exist for only a handful of 
models. Fitting forms, such as in |^, are tuned to the 
cosmological constant model and even then agree with 
simulations to only ^ 10%. 

To use the next generation of large scale structure ob- 
servations we need a better, more efficient method of pre- 
dicting the effects of cosmology on the power spectrum. 
Section ^] presents a new prescription for this, reducing 
the dimensionality of the grid of cosmologies that need 
simulation. In ^IIII we present the results of simulations, 
obtaining ~ 1% accuracy in the matter power spectra. 
The method is extended in i jlVl with a discussion of time 
varying equations of state and cosmic microwave back- 
ground anisotropies. 



II. GOING NONLINEAR 

In general relativity the power spectrum of linear mat- 
ter density perturbations can be readily calculated (see 
4] for a recent convergence study). On large scales and at 
late times the power spectrum grows with fixed shape and 
the growth function can be evaluated, even for dynam- 
ical dark energy cosmologies, from a simple differential 
equation |lj] or highly accurate fitting formulas j5|| . This 
remains true even for some extended theories of gravity 
Since nonlinear structures develop from linear per- 
turbations, the linear growth function is a natural place 
to start. 

Given two cosmological models, we can match their 
growth factors at the present by normalizing to the same 
mass fluctuation amplitude, erg. For definiteness, con- 
sider matching some dark energy model, characterized by 
its equation of state or pressure to energy density ratio, 
w = p/ p, to a, cosmological constant (w — —1) cosmol- 
ogy with dimensionless matter density Qm = 0.3 (we as- 
sume spatial flatness for all models). Next we match the 
normalized growth factors g{a) at a second scale factor. 
This scale factor, a*, is chosen to minimize the deviation 
between the models' linear growth factors over some red- 
shift range. For the range z = — 3, we find a^, = 0.35, or 

— 1.86, is a good choice. One could choose a different 
range and obtain a different a^,. 

Matching the g(a*) determines a value of ri,„ for the 
dark energy model. Finally, we fix the Hubble constant 
by keeping ilmh^ the same for the two models. Table |2 
shows the results for various dark energy cosmologies. 

The degeneracy relation upon changing w, or flm, be- 
tween two constant equation of state models is 

Ar2„ w 0.42(1 - a^)Au;. (1) 

For example, if the models differ by 0.02 in ilm, they 
must also differ by 0.096 in w to match growths at, say, 
z = 1 (generalizing a*). This is only a rule of thumb 
(e.g. a coefficient 0.39 works better for w < —1 models 
and 0.44 for u; > — 1 models), good near the fiducial 
rim = 0.3, w — —1. However it is easy to compute the 
matching numerically, and we use the exact calculation 
throughout. 
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TABLE I: Cosmological parameters generated from our pre- 
scription for matching the nonhnear power spectra of dark 
energy cosmologies. The last column is defined as the present 
day value of the dimensionless power spectrum &t k = H , 
i.e. A^(A:) = 5li{k/ Hf+"T^(k), where we adopt scale invari- 
ance n = 1, and T{k) is the transfer function @, fixed for k 
measured in 1/Mpc. 



III. RESULTS FROM SIMULATIONS 
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FIG. 1; The nonlinear mass power spectrum is shown at three 
redshifts. The top panels plot A^(A;) = k^Pk/{27v') from the 
simulations for five different dark energy models matched ac- 
cording to our prescription. The solid line indicates the Smith 
et al. 1^ fitting formula. The bottom panels show the ratio of 
A^(fc) for the four dark energy models relative to the ACDM 
case. 



By "working at constant flmh^ and flbh^ wc fix the 
shape of the matter power spectrum in physical units, 
e.g. k in Mpc~^. Our procedure arranges a close match 
of the growth of the power spectrum in the linear regime. 
Thus it is natural to expect that the nonlinear power 
spectra will be closely matched^. To test this we make 
use of N-body simulations. 

For the first 5 models of Table HI we ran N-body simu- 
lations using a parallel particle- mesh (PM) code 9]. In 
each case we evolved 256'^ equal mass particles using a 
512'^ force mesh from z = 60 to the present, dumping the 
phase space data at z ~ 1.86, 1, and 0. The simulation 
volume for the results presented here was a periodic, cubi- 
cal box of side 366 Mpc, though we also ran a few larger 
boxes to check convergence on the largest scales. For 
each output time, the power spectrum was computed by 
assigning the particles to the nearest grid point of a reg- 
ular, 512^ Cartesian mesh and Fourier transforming the 
resultant density field. The resulting A^(fc), corrected 
for the assignment to the grid using the appropriate win- 
dow function and for Poisson shot noise, were placed in 
logarithmically spaced bins of k. The average A^(fc) is 
plotted at the position of the average k in each bin. 

While the PM code has limited force resolution, it al- 
lows high mass resolution at modest computational cost. 
This is important in the translinear scales (A^ 1) of 
interest here. The finite force and mass resolution of the 
PM code limits the dynamic range in fc to ~ 50, and the 



short-fall in power can be easily seen in the highest k 
point plotted in the top left panel of Fig.^ Much of the 
effect of the force resolution is reduced in taking the ratio 
(see also [l^ ) of power spectra, as in the lower panels of 
Fig. n although the details of the shot-noise correction 
are starting to become important on these scales. 

From Fig. ^ we see that across most of the range of 
scales probed by, e.g., weak lensing^ and for which we 
believe purely gravitational N-body simulations will pro- 
vide accurate calculations of the power spectrum, the 
models are indeed nearly degenerate. The maximum de- 
viation in A'^{k) is 1.5%, for the fim = 0.26, w = —1.156 
model at z = 0.5. This is also the maximum deviation 
in the (square of the) linear theory growth rates. Where 
the growth rates are closely matched, the nonlinear power 
spectra are too. 

IV. GROWTH AND THE CMB 

A significant new result is that the distance to last scat- 
tering diss is automatically preserved through our match- 
ing prescription. That is, by matching the growth today 
and at a* = 0.35, we ensure that diss changes very lit- 
tle; the growth factor determines the distance. The frac- 
tional distance deviation of the dark energy models from 



^ This is certainly true to second order in perturbation theory 
— see lal for a review. 



^ For sources at Zs ~ 1, or comoving distance 3Gpc, wavenum- 
ber k = lMpc~^ corresponds to multipole £ ~ kx = 1500, or 
arcminute scales, at the peak of the lensing kernel, x = 1.5 Gpc. 
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the ACDM case is less than 0.04% (e.g. 0.001% for the 
= 0.28, w = -1.076 case). 

This is important because the CMB is one of our 
strongest probes of large scale structure and cosmology. 
The main impact of the dark energy and the matter den- 
sity on the CMB power spectrum is through the geomet- 
ric degeneracy involving the distance to the last scatter- 
ing surface, and the physical matter density Vlmh? . Our 
prescription takes care of both, meaning that we select 
models that automatically include a strong CMB prior 
(see also 0): our matched models, with O^/i^ and diss 
fixed, will agree on the CMB temperature power spec- 
trum to high accuracy, except at low multipoles where 
the integrated Sachs- Wolfe effect enters (though cosmic 
variance makes differences hard to detect) and at high 
multipoles through the effect of gravitational lensing. 

We can further extend this to models with time varying 
dark energy equation of state. For values of wq and Wa, 
where w{a) = wo + Wa(l — a), we match the growth at the 
present and at a^, (we have checked that a* = 0.35 works 
for these cases as well) . This provides the necessary value 
of Om, and hence h (see Tabled). For matching models 
with {wo.Wa) of (-1,0.3), (-1,-0.3), and (-0.8,-0.3), 
the maximum deviation from the fiducial ACDM model 
in the linear growth factor is less than 0.6%, dX z — 0.6. 
So we expect, based on illlll that the nonlinear power 
spectra will match to within 1.2%. A matching model 
with (—0.8, 0.3), where the average equation of state lies 
further from the fiducial w = —1 and the required 
is 30% different from the fiducial model, fits somewhat 



less well: the expected deviation in the nonlinear power 
spectrum is 3.4%. 

The time varying equation of state models still preserve 
the CMB matching as well. They match on diss to better 
than 0.1%, except for (—0.8,0.3) which agrees to 0.4%. 



V. CONCLUSION 

We present a simple prescription for generating the 
nonlinear mass power spectrum of dark energy cosmolo- 
gies, accurate to ~ 1 — 2%. It requires only straight- 
forward calculation of the linear growth factor at two 
redshifts. This will allow efficient generation of a suite of 
dark energy cosmologies using only a reduced dimension 
grid of ACDM simulations. 

Moreover, we have identified a strong degeneracy be- 
tween the growth factor and the CMB distance to last 
scattering, such that our matching procedure can auto- 
matically satisfy CMB constraints. 
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